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A CO (3-2) survey of a merging sequence of luminous 
infrared galaxies 
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ABSTRACT 



Luminous infrared galaxies (£tr, > 10 £©) are often associated with interact- 
ing galactic systems and are thought to be powered by merger-induced starbursts 
and/or dust-enshrouded AGN. In such systems, the evolution of the dense, star form- 
ing molecular g function of merger separation is of particular interest. Here, 
we present observations of the CO (3-2) emission from a sample of luminous infrared 
galaxy mergers that span a range of galaxy-galaxy separations. The excitation of 
the molecular gas is studied by examining the CO(3-2)/CO(l-0) line ratio, r^i, as 
a function of merger extent. We find these line ratios, rax, to be consistent with ki- 
netic temperatures of 7i-=(30-50) K and gas densities of hh 2 = 10 3 cm -3 . We also 
find weak correlations between rai and both merger progression and star formation 
efficiency (LfiR,/Lco(i-o))- These correlations show a tendency for gas excitation to 
increase as the merger progresses and the star formation efficiency rises. To conclude, 
we calculate the contributions of the CO(3-2) line to the 850 /im fluxes measured with 
SCUBA, which are seen to be significant (~24%). 

Key words: galaxies: starburst - galaxies: interactions - galaxies: ISM - infrared: 
galaxies - ISM: molecules 



1 INTRODUCTION 

Luminous Infrared Galaxies (LIGs) (Lir > 10 11 Lq) and Ul- 
tra Luminous Infrared Galaxies (ULIGs) (Lir > 1O 12 L0) 
are the dominant population for galaxies with bolomet- 
ric luminosi ties greater than 1 1 Lq in the local universe 
(z < 0.3) (|Soifer et al.l Il986l h Observations have shown 
that most are associated with mergin g/interacting galax- 
ies wh ich are rich in molecular gas fsee lSanders fc Mirabell 
(1996) for an extensive review). Analysis of the original 
Infrared Astronomical Satellite (IRAS) data showed that 
LIGs are quite rare in the local u niverse but increase 
significantly with inc reasing redshift ( Hac king et alj 1 19871; 
IChapman etaH l2005l h Studies of local luminous infrared 
galaxies have, therefore, become essential for a wider under- 
standing of the role of mergers and starbursts both in the 
local universe and at higher redshifts, where such objects 
are more commonly found. 

Optical studies of local LIGs have shown that many 
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are associated with interacting or merging galactic systems 



llArmus et al 



19871: IScoville et al.ll 19981: lArribas et al]|2000l : 



IVeilleux et al. 2002j). It is now thought that as a galactic 
merger progresses, gas and dust from the parent galaxies 
dissipate energy through shocks, giving rise to an in-fall of 
gas and dust towards the centre of gravity of the interacting 
system l|Barnes fc Hernquistll 996 ; Mih os fc Hernquistll996l : 



ICox et al.ll20Q6h . The concentration of cold, dense material 
then acts to fuel starburst activ ity, or possibly an active 
galactic nucleus (AGN), or both i|Kim et alj|l998h . The en- 
ergy produced by the starburst or AGN in the optical and 
ultraviolet is absorbed and then re-radiated by dust, leading 
to an intense infrared luminosity. 

The effect of a galactic merger upon star formation, 
as well as its effect upon the total mass and properties of 
the molecular and atomic gas within the interacting sys- 
tem , has been the s ubjec t of several observational stud- 
ies. iGao fc Solomon! l| 19991 ) used the CO (1-0) spectral line 
strength as an indicator of total molecular gas mass, Mh 2 
and plotted this quantity against the projected distance of 
the merging galactic nuclei at optical wavelengths, the pro- 
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jected component separation, in a merging sample of LIGs. 
The authors found a positive correlation between these two 
quantities. The authors also found an anti-correlation be- 
tween the ratio of infrared luminosity to molecular gas mass, 
Lir/JV/h 2 , a measure of star formation efficiency, and the 
component galaxy separation for their sample. Both obser- 
vations are consistent with a depletion of gas as a huge 
increase in star formation (to ar ound ~ 100 Mp, yr~ ) oc - 
curs. In a survey of literature data. lGeorgakakis et al.l (|2000l ) 
found strong evidence for gas depletion and increased star 
formation efficiency in a merging galaxy sample which in- 
cluded pre-merger and post-merger candidates. The au- 
thors also found an increase in the central molecular hydro- 
gen surface density, as traced by CO (1-0), and a decrease 
in the fraction of cold gas mass as the galactic mergers pro- 
gressed. This result is consistent with large molecular gas 
inflows, as predicted by numerical simulations, the conver- 
sion of neutral hydrogen to molecular hydrogen, stars and 
hot gas and molecular hydrogen depletion due to ongoing 
star formation. 

Merging and starburs t activity also affec ts the atomic 
component of the ISM. iMirabel fc Sanders! |l989) have 
found that the ratio of molecular hydrogen to atomic hy- 
drogen gas mass, Mh 2 /Mhi, increases with Lir. It is cur- 
rently unclear whether this increase is chiefly due to a de- 
pletion of HI via ejection from the interacting system, de- 
pletion of HI due to photo-ionisation caused by the for- 
mation of young stars or enhanced formation of molecular 
clouds in merger induced shocks. Further observations of 
Mhi have been attempted to investigate whether this deple- 
tion of atomic gas is also se en across a merging sequence of 
LIGs (Ivan Driel et al.ll200ll ). but no conclusive answer has 
yet been derived. 

The evolution of molecular hydrogen gas, which is as- 
sociated with star formation, within merging luminous in- 
frared galaxies is of particular interest. Observations of the 
CO(l-O) rotational transition are often used to trace the 
total molecular hydrogen gas mass e.g. lYoung fc Scovilie] 
(1991). After H2, CO is the most abundant molecule in 
the ISM ([CO/H2] ~ 10" 4 ), and the CO(l-O) transition 
is typically thermalised in molecular clouds, making it a 
very convenient transition to observe. The low excitation en- 
ergy, Eio/ks = 5.5 K and the low critical density 7iH 2 ,crit ~ 
410 cm -3 of the CO(l-O) transition make it ideal for tracing 
the bulk metal-enriched H2 in galaxies. These same prop- 
erties, however, make CO (1-0) insensitive to the physical 
conditions of the molecular gas. If we require constraints 
on the density and kinetic temperature, we can o bserve the 
highe r- J CO lines, (i.e. COj+i-j, where J ^ 2) (|Zhu et al.l 
120031 ). These transitions will only become thermalised and 
luminous for the denser and warmer gas components and it 
is these dense components which are a direct indication of 
the star formation potential of the interacting systems. 

The higher order CO (3-2) transition, in particular, is 
commonly used to trace the warmer, d enser components of 
the I S M associated with star formation (IMauersberger et aU 
Il999l : iBavet et all l200d : IWilson et aU l2009h and constrain 
the kinetic temperature and density of the molecular gas. 
It is likely that there will be multiple density compo- 
nents of molecular gas within the ISM of luminous infrared 
galaxies, and t he critical de nsity of CO (3-2) (nH 2 , cr it > 
8.4 x 10 3 cm" 3 . Hansen! I1995T I is well matched to the den- 



sity of the star forming component. The relative ratios 
of the higher-order CO line intensities depend on both 
the density and the kinetic temperature of the molec ular 
gas jScoville fc Solomon|[T973 ; lGoldreich fc K wad! 1974 ). By 
measuring a range of CO line ratios one can therefore begin 
to constrain the likely conditions of the molecular interstel- 
lar medium. Observing the CO (3-2) transition is convenient 
as it traces mostly the starforming molecular gas and has 
relatively large luminosities compared to other commonly 
used tracers, such as HCN (nH 2 ,crit > 10 s cm -3 ), which have 
higher dipole moments, but lower abundances in the ISM. 
Measuring the CO (3-2) emission is also important since 
this emission can contaminate 850 /im continuum observa- 
tions of the dust component on the ISM for these galaxies 
l|Papadopoulos fc Allen! [2000 ) . Correcting this 850 /im flux 
for CO emission is important when using thi s flux to deter- 
mine the dust mass and kinetic temperature (jSeaquist et al.l 
l2004h . 



In this paper, we present observations of the CO(3-2) 
rotational transition to investigate the evolution of the dense 
(^H 2 ,crit > 8.4 x 10 3 cm -3 ) molecular gas component as a 
function of merger extent in a sample of LIG mergers. Our 
aim is to study the evolution of merging luminous infrared 
galaxies by observing a LIG sample at differing stages of in- 
teraction, and thus obtain a detailed understanding of how 
the starburst evolves as the interaction proceeds. Such obser- 
vations are key to obtaining a temporal understanding of the 
evolution of the molecular gas reservoir and how it relates to 
star formation activity as the merging component galaxies 
become increasin gly tidally disrup t ed. O ur study builds on 
work reported bv lGao fc Solomon! (|l999T ). in which a sam- 
ple of luminous galaxies was chosen to represent a merging 
sequence with a range of nuclear separations of the merging 
components. We observe the CO(3-2) transition in a subset 
of this merging sample of 49 LIGs to trace the evolution 
of the warmer, denser molecular gas component as a func- 
tion of merger separation. In particular, we investigate how 
the excitation of the molecular gas changes by examining the 
evolution of the CO (3-2) /CO (1-0) line ratio, r^i, as merging 
progresses. r$i will be larger for warmer, denser gas and one 
might expect its value to increase with decreasing nuclear 
separation, as an increasing fraction of the total molecular 
gas becomes concentrated in star forming regions. We then 
examine how the line ratio varies with merging extent, far- 
infrared luminosity (Lar, defined in Section [2]) and star for- 
mation efficiency (Lfm/Lcon— o))> an d compare our sample 
with two samples with smaller infrared luminosities. Finally 
we determine contamination arising from the CO (3-2) line 
flux upon the 850 ^m continuum observations of the dust 
component of the ISM for these galaxies. Section [2] outlines 
the observations that were made, including the sample selec- 
tion criterion, the choice of pointings and calibration of the 
data. The data reduction and results of the observations are 
presented in Section [3] and their interpretation and correla- 
tions with other data for the sample are discussed in Section 
[4] Section [5] outlines the conclusions which may be drawn 
from this study. 
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2 SAMPLE SELECTION AND OBSERVATIONS 

2.1 The sample 

Our sample of merging gal axies is a subset of a s ample of 
merging galaxies defined bv lGao fc Solomonl j 19991 ). chosen 
such that the sources are ob servable by the James Clerk 
Maxwell Telescope (JCMTlFl lGao fc Solomonl (|l999h chose 
a parent sample of 49 objects to include all LIGs with (a) 
available CO (1-0) and R-band CCD data (b) i m ) 2 x 
lO 11 !/^, and (c) having nuclear separations (S scp ) greater 
than 2" and less than half of the sum of the galactic op- 
tical disc sizes. The constraint on nuclear separations was 
set to ensure that all genuinely interacting systems were in- 
cluded in the sample, whilst excluding difficult to measure, 
late mergers with separations of S scp < 2", as they are in- 
distinguishable from merged single-nucleusgalaxies with the 
usual seeing limit for optical observations □ 

The sample chosen for this study, hereafter referred to 
as the CO (3-2) sample, comprises of 33 objects from the Gao 
& Solomon sample, chosen to meet the following selection 
criteria: 

(i) A declination such as to be observable from the James 
Clerk Maxwell Telescope (S > -40°). 

(ii) A distance such that the redshifted CO (3-2) emis- 
sion line falls within the 315-373 GHz tuning range of the 
345 GHz receiver, and at least 2 GHz either side of the strong 
325 GHz atmospheric water vapour absorption feature (i.e. 

< cz < 17230 km/s and 21160 < cz < 29310 km/s). 

Shown in Figure [1] is the statistical spread of nuclear 
separation, infrared luminosity and far-infrared luminosity, 
for the Gao & Solomon sample and the CO (3-2) sample 
defined here. 

2.2 Choice of pointings 

The beam width (FWHM) for Receiver B3 on the JCMT is 
14" and the nuclear separations for our sample lie between 
2.5" and 40.6", which made it necessary to choose suitable 
pointings carefully so as to capture all of the CO (3-2) emis- 
sion. The choices of pointing centre were made based on 
nuclear separation, which in turn were deri ved from either 
images taken at R-band from I Sanders! ll 19921 '). 645 nm (band 
103aE) ima ges from the STScI Digitized Sky Survey (DSS) 
(|Kentlll994h or with WFPC2 images taken with the Hubble 
Space Telescope. Where available, high-resolution CO (1-0) 

1 The JCMT is run by the UK Science and Technology Facil- 
ities Council, the National Research Council, Canada, and the 
Netherlands Organisation for Scientific Research. 

2 Infrared luminosities follow standard IRAS definitions: Li R = 
L(8 - lOOO^m) = 47rD 2 j F IR [L ] and L an = L(40 - 500Atm) = 
47rD 2 CF fIR [L Q ] where F m = 1.8 x 10- 14 (13.48/i 2 + 5.16/ 25 + 
2.58/eo + /i oo)Wm- 2 and F flR = 1.26 x lQ- 14 (2.58/ 60 + 
flop) Wm~ 2 JSanders fc Mirabej|l99l|PerauhJll987l ; lHelou et al.l 
11988) ■ 

3 Since the publication of iGao fc Solomoiil il999h , the nuclear 
separation measurement for IRAS 17208-0014 has been improved 
and the object has been found to be a late merger (Sscp ^ 2 ) 
we thus exclude it from all further sub-samples discussed in this 
paper. 



interferomet ry maps were also used to confirm the choice 
of po intings (|Wang et al.ll200ll : ICasoli et al.|[l999al : iLo et ail 
1997). Three different categories of pointings/observations 
were identified: 

(i) Objects with a projected nuclear separation of < 7" 
(less than half of beam size): single pointing sufficient to 
capture all of the CO(3-2) emission. 16 objects of our sample 
were observed with a single pointing. 

(ii) Objects with a projected nuclear separation of ~ 14" 
(of order of the beam size): two pointings used, centred on 
the nuclei of the merging galaxies. 7 objects of our sample 
were observed with a dual pointing. 

(iii) Objects with a projected nuclear separation > 14": 
multiple pointings, separated by 1/2 beam (7") were used to 
realise a fully sampled strip along the axis joining the two 
component nuclei. 5 objects of our sample were observed 
with multiple pointings (see Figure [CT]) . 

2.3 CO(3-2) observations 

Observations of the CO(3-2) line (y reBt = 345.796 GHz) 
were made using the dual-channel (orthogonal polarisation), 
345 GHz receiver (RxB3) on the JCMT. The DAS (Digital 
Auto-correlation Spectrometer) backend was used in one of 
two bandwidth modes, depending on the expected width of 
the CO(3-2) spectral line given the CO(1-0) linewidths. The 
single-channel, wideband mode, in which the two 900 MHz 
subbands are concatenated to realise a useable contiguous 
bandwidth of 1800 MHz (Av ~ 1560 km/s), was used for 
expected line widths of Afwhm > 500km/s; the more sen- 
sitive dual-channel standard mode, with a single subband 
of 900 MHz (Av ~ 780 km/s) was used where lines of 
Awpwhm < 500km/s were anticipated. The wideband mode 
was only used where the velocity covered by standard mode 
was not sufficient to provide line-free regions for baseline 
removal. 

Rapid beam-switching at a secondary mirror chop fre- 
quency of 1-2 Hz and beam throws of 60 — 90" in azimuth 
were used to ensure good atmospheric cancellation. Hot/cold 
load calibrations were performed every 20 to 30 minutes to 
achieve consistent calibration to the T* A antenna tempera- 
ture scale. The telescope pointing was checked every 2-3 
hours using local spectral line pointing sources, and any off- 
set corrections applied to the telescope pointing model, re- 
sulting in an RMS pointing accuracy of 1-2". The science 
data taken for each pointing consisted of spectra from each 
RxB3 polarisation channel, usually integrated for a period of 
around 10 minutes. Integration times of longer than 10 min- 
utes were achieved by coadding several ~ 10 minute spectra 
while checking the data for instrumental drifts and baseline 
problems. 

Spectral line standards such as OMC1 were observed to 
check amplitude calibration and quantify calibration uncer- 
tainties, as well as to monitor the overall performance of the 
telescope. Mars and Uranus, when observable, were observed 
to establish main-beam efficiencies. 

A total of 28 of the 33 sources in the CO (3-2) sample 
were observed over the period between October, 2000 - July, 
2003. The atmospheric opacity at 225 GHz measured at the 
CSO, was typically T225GHz ~ 0.1 giving system tempera- 
tures of around 600 K in the receiver tuning range used for 
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Figure 1. Nuclear separation (kpc), infrared luminosity Ljr and far infrared luminosity Lfia statistics for the galaxies in the Gao and 
Solomon sample (open rectangles), the fraction of the sample observable at the JCMT (narrow cross-hatch) and the 28 galaxies observed 
in this study (wide cross-hatch). 



the observations. Of our 28 observed sources ,16/7/5 
were observed with single / two / multiple pointings. Fig- 
ure [T] shows a comparison of the statistical spread of nuclear 
separations, infrared luminosities and far-infrared luminosi- 
ties for the objects actually observed. 



3 DATA REDUCTION AND RESULTS 

All data were reduced u sing the JCMT spectral line re- 
duction package SPECX |PadmaTj|l993h . Individual spec- 
tra were inspected, with any scans showing significant base- 
line structure/ripple rejected (~ 5% of total). The spectra 
were then coadded, linear baselines fitted and removed and 
then the spectra were binned to 26 km/s (30 MHz). The 
Ico(3 — 2) values were determined by directly integrating 
under the extent of the spectral line. 

Planetary observations of Mars and Uranus were used 
to estimate »7 m b and calibrate the data directly from the 
sky-corrected antenna temperature, T% to the main-beam 
brightness temperature T m b via 

7? mb = B eB /Fes = ^r~. (1) 

J-mh 

Here, -B e ff, the beam efficiency, and _F e ff, the forw ard effi- 
ciency follow the definitions given in lDownesl fl98^). Values 
of »7 m b ~ 0.55 were typically measured using this technique. 
Where no planetary observations were available, a canonical 
value for RxB3 on the JCMT of ?7 m b = 0.55 was adopted, 
which reflects the mean value obtained and published on the 
JCMT web pages for the observing perioqj. The observed 
spread for 7/ m b over the period of observing was ~ 20%, 
consistent with that typically seen for RxB3 on the JCMT. 

The uncertainty in 7/ m b introduces a systematic calibra- 
tion uncertainty of ± 20%. In addition there are random 
noise contributions arising from (a) noise across each spec- 
tral channel (b) determining the integrated spectral line and 
(c) fitting the overall subtr acted basel ine level. These errors 
can be combined, following lGad (| 1996T) (see AppendixfS]), to 
give an overall error, <r/co(3 — 2), in the velocity-integrated 
main-beam brightness temperature (the integrated line in- 
tensity), 7co(3 — 2). 



: //www. jach .hawaii . edu/ JCMT/ spectral. line/St andards/beameff .html 



The efficiency-corrected spectra together with the 
pointing centres superimposed on a Digital Sky Survey op- 
tical image, are shown in Figure [ClI at the end of the paper. 
The properties of the spectra for each pointing, together 
with the random errors, ol are summarised in Table [1] 
Where no line was detected, an upper limit for 7co (3 — 2) of 
2cr is given, assuming the detected line would have had a full 
width to zero intensity of AW = 400 km/s. The observed 
spectra shown in Figure IC1I show a wide range of effective 
linewidths (from 74 to to 639 km/s) and often show double 
peaked velocity profiles. These profiles can arise when a sin- 
gle telescope beam picks up CO (3-2) emission from each of 
the interacting nuclear components or when there is rotation 
of the molecular gas around one of the nuclear components 
within the beam. 

3.1 Luminosities and line ratios 

Having determined 7co(3 — 2) for each pointing for each 
object, we proceeded to calculate the total CO (3-2) line lu- 
minosity and the CO(3-2)/CO(l-0) line luminosity ratio for 
each object as follows. For galaxies with a projected nuclear 
separation of < 7", it was assumed that the angular size 
of the CO (3-2) emitting region was small compared to the 
14" beam size w hen determining the ove rall CO (3-2) lumi- 
nosity. Following [Gao fc Solomon] (|l999l ). the CO(3-2) line 
luminosities, 7/co(3-2), were calculate d by following the def i- 
nition of CO (1-0) luminosity given bv lSolomon et all i|l997h . 

^CO(3-2) = TbAVD. s D\ 

= 23.5fU b Di/ CO (3-2)(l + ^r 3 [Kkms-'pc 2 ] (2) 

where Tb is the peak line brightness temperature, AV is the 
effective line width, 7qo(3-2) = / TmhdV, T m b is the main 
beam brightness temperature, Dl is the luminosity distance 
and is the solid angle (arcsec 2 ) of the source convolved 
with the beam. A value of Q. sir b = fib = 1.133t9f whm — 
222 arcsec 2 was used for objects observed with a single point- 
ing, following the above assumption that the CO (3-2) emit- 
ting region was small compared to the beam. For objects 
observed with multiple pointings the same value was used 
after combining the individual fluxes from each pointing, 
^co(3-2), to calculate a total flux, ico(3-2,Totai)- Upon ex- 
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Table 1. Observational parameters for the LIGs observed in this study. 
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(0,8) 




1.7 
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(-!,-') 




2.6 












(-2,-14) 


< 5.2 


2. 












(-3,-21) 


< 4.4 


2-2 












(-3,-28) 


< 5.2 


2.6 












(-4,-37) 


8.79 




74 


Mrk 848 


15:18:06.2 


+42:44:42 


12043 


(0,0) 


19.9 


1.1 


121 


NGC 6090 


16:11:40.7 


+ 52:27:25 


8754 


(0,0) 


56.4 


2.2 


141 


17132+5313 


17 14 20.38 


+ 53 10 30.4 


15290 


(0,0) 


12.5 


2.5 


221 


NGC 6670 


18:33:37.7 


+ 59:53:24 


8684 


(0,0) 


13.9 


1.2 


151 










(-41,-15) 


<3.8 


1.9 












(-35,-11) 


11.4 


1.9 


166 










(-28,-9) 


<4.0 


2.0 












(-21,-6) 


19.2 


2.5 


245 










(-14,-3) 


12.1 


1.8 


174 










(-7,-2) 


5.7 


1.0 


100 


19297-0406 


19:32:21.3 


-03:59:54 


25674 


(0,0) 


<1.8 


0.9 




20010-2352 


20:04:03.5 


-23:43:54 


15249 


(0,0) 


16.4 


1.4 


344 


II Zw 96 


20:57:23.7 


+ 17:07:44 


10900 


(0,0) 


13.7 


1.4 


176 










(8,-9) 


17.9 


1.4 


180 


NGC 7592 


23:18:22.6 


-04:24:58 


7328 


(0,0) 


<5.2 


2.6 












(-12,2) 


<6.2 


3.1 




Optical rcdshift definition 


, values 


obtained 


from 


NASA/IPAC 


Extragalactic 


Database (NED) 



|http: //n edwww . ipac . caltech . edu/ | 
Effective linewidth = T^b.peak /IcO (3 — 2), where T^t^peak is the peak main— beam brightness temperature of the line. 



amining SCUBA 850/im continuum maps for Mrk 1027 after 
we made our CO (3-2) observations we found that the dust 
emission from the source appeared extended (6" by 12") 
and that centre of this emission was offset by 3.5" relative 
to our CO(3-2) (0,0) pointing. We thus took /co (3 -2, Total) = 
1.6* 1.4/(0, 0), with the first factor correcting for beam dilu- 
tion and the second factor correcting for the pointing offset. 
For IRAS 10039-3338 the flux was inadvertedly measured 
at 3.5" from centre of optical/ IR emission, so we corrected 



this mispointing by using 7 C o(3-2,Totai) = 1.47(0,0), which 
assumes the emission is small compared to the beam. 

For galaxies observed by more than one pointing it was 
necessary to combine the fluxes to determine the total line 
luminosity. For galaxies with a projected nuclear separation 
~ 14" observed with two pointings, the fluxes were com- 
bined as a weighted sum using the method presented in Ap- 
pendix [B] For objects with nuclear separation greater than 
14", the fluxes from each pointing were combined by assum- 
ing a plausible underlying brightness distribution based on 
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Table 2. Component separations, infrared and CO luminosities and line ratios for our LIG sample. 



Object 



Separation 



ico(3-2) 



Tico(3-2) L C o(3-2) Loo(l-O) 





(kpc) 


(") 


1QllL e 


1q11l 


(K km/s, Total) 


(K km/s) 


(10 9 L,) a 


(10 9 L,) a 




00057+4021 


2.1 


2.5 


4.3 


3.31 


9.5 


0.5 


1.65 


3.8 


0.43 


Arp 236 


(i.l 


16 


4.71 


3.19 


152 


3.0 


5.21 


10.67 


0.49 


01077-1707 


22.4 


35 


4.33 


2.79 


24.6 


2.4 


2.49 


8.5 


0.29 


III Zw 35 


4.1 


8.2 


4.15 


3.58 


18.5 


2.2 


1.21 


2.09 


0.58 


Mrk 1027 


5 


9 


2.49 


1.69 


63.8 


2.6 


4.79 


5.35 


0.89 


02483+4302 


3.6 


3.8 


6.21 


4.54 


5.74 


0.7 


1.30 


2.88 


0.45 


UGC 2369 


13.1 


22.5 


3.95 


2.66 


20.5 


2.3 


1.69 


7.25 


0.23 


03359+1523 


6.5 


10 


3.34 


2.42 


12.3 


2.5 


1.30 


6.9 


0.19 


04232+1436 


6.3 


4.6 


11.15 


7.49 


20.0 


2.4 


10.3 


9.0 


1.15 


Arp 55 


10.7 


15 


4.66 


3.39 


45.1 


3.8 


5.74 


12.07 


0.48 


10039-3338 


3.1 


5 


4.54 


3.94 


15.3 


0.6 


1.46 


2.93 


0.50 


10190+1322 


6.6 


5 


11.34 


7.13 


10.3 


1.1 


4.86 


8.2 


0.59 


10DOO + 244ft 


6.2 


8 


9.58 


7.59 


28.4 


0.0 


4.31 


5.76 


0. 75 


Arp 299 


4.5 


25 


6.36 


4.74 


183 


2.8 


1.64 


3.0 


0.55 


13001-2339 


2.6 


6.5 


2.45 


2.0 


43.8 


1.5 


1.67 


2.55 


0.66 


Arp 238 


12.1 


21 


5.17 


3.8 


16.3 


2.7 


1.34 


4.24 


0.32 


NGC 5256 


4.8 


9.5 


3.05 


2.02 


39.6 


4.2 


2.46 


5.66 


0.43 


Mrk 673 


4.0 


6 


2.2 


1.51 


18.9 


2.8 


2.04 


4.0 


0.51 


14348-1447 


5.7 


4.0 


20.43 


16.28 


19.4 


2.0 


10.55 


14.0 


0.75 


Arp 302 


25.8 


41.6 


4.1 


3.59 


59.0 


2.4 


5.59 


19.7 


0.28 


Mrk 848 


4.8 


6.5 


7.16 


5.29 


19.9 


1.1 


2.62 


7.03 


0.37 


NGC 6090 


3.5 


6.5 


2.98 


1.94 


56.4 


2.2 


3.96 


5.0 


0.79 


17132+5313 


6.4 


6.5 


7.67 


5.73 


12.5 


2.5 


2.64 


7.24 


0.37 


NGC 6670 


14.6 


26.5 


3.84 


2.73 


44.6 


3.4 


3.26 


11.6 


0.28 


19297-0406 


11.5 


8 


24.59 


19.28 


<1.8 


1.1 


<1.1 


9.43 


<0.11 


20010-2352 


8.1 


8.9 


4.69 


3.16 


16.4 


1.4 


3.47 


6.38 


0.54 


II Zw 96 


7.4 


11 


7.57 


6.36 


27.2 


1.7 


3.08 


6.04 


0.51 


NGC 7592 


5.3 


11.5 


2.43 


1.68 


<6.2 


3.1 


<0.31 


4.55 


<0.07 


Mean 


7.15 


11.5 


5.2 


4.03 


37.8 


2.0 


3.2 


6.5 


0.47 


Calculated using equation 


11 , H 


— 75 km/s 


Mpc" 1 , 90 


=0.5, L t = Kkms" 


-!pc 2 . 









the object morphology. For IRAS 01077-1707 the emission 
appeared consistent with a single compact source at (14, 0), 
so we estimated the total flux via 7co(3-2, Total) = 7(14,0). 
For Arp 236, the emission of the source appeared to be 
consistent with a source equal to the beamsize so we took 
^CO(3-2, Total) = \/2~Zco(3 - 2)(0,0). For Arp 299 we com- 
bined the fluxes via 7 C 0(3-2,Totai) = (1(0,0) + I(-19,4) + 
J(— 7, 1) + /( — 14,2))/2, consist ent with an i nterfer omet- 
ric CO(1-0) map appearing in ICasoli et al] l|l999bf ). For 
Arp 302 we assumed a uniform underlying distribution 
across the pointings at (0,0), (0,8) and (-1,-7) and combined 
these separately with the flux at (-4,-37), /co(3-2, Total) = 
(1(0, 0)+/(0, 8)+I(-l, -7))/2+I(-4, -37). For NGC 6670, 
the emission appeared to be strongest at pointings (0,0), 
(-21,-6) and (-35,-11), so we estimated the total flux via 
/co ( 3-2,Totai) = 1(0, 0) 4- J(-21, -6) 4- J(-35, -11), i.e. as- 
suming that most of the emission is concentrated in these 
directions. The derived values for the total CO(3-2) lumi- 
nosity ico(3-2) an d the CO(3-2)/CO(l-0) line luminosity 
ratios, r3i are presented in Table [2] 

The full physical implications of line ratio measure- 
ments for molecular gas conditions require non-LTE pho- 
ton transport models. These can be solved numerically, ac- 
counting for the excitation between several energy levels of 
the molecular species present. A popu lar approach is the 
Large Velocity Gradient method (LVG) IjScoville fe Solomon! 
1 19741 ; iGoldreich fc Kwanl [l97i l which assumes that large- 
scale systematic velocity gradients, rather than thermal mo- 
tions, dominate the observed line widths. This assumption 
simplifies the photon transport problem, since there is es- 
sentially no overlap in line emission between distant parts 
of the cloud. For the mean value of rsi — 0.47 for the 



LIGs in our sample, we performed an LVG fit assuming 
a 12 CO(1-0)/ 13 CO(1-0) line ratio, R 10 of .Bio = 15-20 
measured typically for such gala xies fe.g. lCasoli et al]|l992l ; 
iPapadopoulos fc Seaquistlll998l ). This LVG fit yielded typ- 
ical conditions of kinetic temperature T k — (30-50) K and 
a H2 number density of n = 10 3 cm~ 3 with a CO abun- 
dance per velocity gradient A = [CO/H 2 ]/(dV/dR) of 10" 5 
(kms -1 pc -1 ) -1 , though good solutions exist also for a much 
warmer (T k ~ 80-110 K) and diffuse (n ~ 3 x 10 2 cm -3 ) gas 
phase. This degeneracy reflects both the lack of constraints 
that could be set by observing a larger number of molecular 
lines (especially high- J CO transitions, 13 CO lines, and high 
density tracers such as HCN transitions), and the presence of 
a diffuse and warm gas phase along with a dense r and cooler 
one t hat usually exists in such systems (e.g. lAalto et al.l 
Il995l ). The hereby observed CO (3-2) line is the highest- J CO 
transition to have significant contributi ons from the former 
while tracing mostly the later phase (Papadop oulos et al.l 
l2007h . We intend to combine our CO(3-2) data with HCN 
observations for the same galaxy sample to further constrain 
the molecular gas conditions. We also intend to conduct a 
more detailed analysis of several of the LIGs in our sample, 
as part of a multi-J CO and HCN line survey of such systems 
l|Papadopoulos et ai]|2007l ). 



3.2 Statistical properties of the sample 

It is instructive to compare the CO (3-2) observations of the 
LIG sample here with the CO (3- 2) observations of a s ample 
of 29 local galaxies reported bv iMauersberger et all (|l999h 
(Figures [2] and [3)l . The authors selected galaxies which ex- 
hibit strong 7co(l — 0) and 7co(2 — 1) intensities, though 
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Number of Galaxies vs L„ 



Number of Galaxies vs r 31 





Figure 2. A comparison of the LfjR for the SLUGS sample (open 
rectangles), our sample (narrow cross-hatch) and the Mauers- 
berger sample (wide cross-hatch). 



Figure 3. A comparison of the rsi for the SLUGS sample (open 
rectangles), our sample (narrow cross-hatch) and the Maucrs- 
berger sample (wide cross-hatch). 



mostly due to their proximity rather than high intrinsic 
CO(l-O) or C0(2-l) luminosity. Their sample contained 27 
non-infrared luminous galaxies with 0.9 x 10 9 Lq < Lib < 
96 x 1O 9 L and 2 LIGs, NGC 2146 (L m = 1.1 x lO 11 !^) 
and Arp 220 (Lm = 1.3 x 1O 12 L ). The majority (26) of the 
galaxies were isolated spirals of type SAB, SA or SB. The 
CO(3-2)/CO(l-0) line ratio for most of the galaxies was be- 
tween 0.2 and 0.7 with only 4 objects exhibiting line ratios 
greater than unity. Thus, most of our high infrared lumi- 
nosity sample have a similar range of CO(3-2)/CO(l-0) line 
ratios to the Mauersberger sample i.e. a selection of nearby 
galaxies with moderate infrared luminosity. 

We can also compare our CO (3-2) observations with 
those of a subsample of the SC UBA Local Univer se Galaxy 
Survey (SLUGS) presented by lYao et all (|2003T ). The 60 
galaxies chosen were a near-complete, flux-limited sub- 
sample Sdo^m > 5.24 Jy sample with Lar predominantly 
less than 10 11 Lq. Their sub-sample is therefore on aver- 
age considerably less IR- luminous than our sample CO (3-2) 
(Figure [2]), yet we see from Figure [3]that a similar spread of 
line ratios is observed with the majority exhibiting < 1. 



4 DISCUSSION 

4.1 Evolution across the merging sequence 

Since our sample represents a merging sequence of galaxies, 
we can now examine how the dense molecular gas compo- 
nent, traced by the CO(3-2) transition, evolves as merging 
progresses. In particular we will determine the correlation 
between the CO(3-2)/CO(l-0) line ratio, m, and nuclear 
component separation, Lan and -i/fiR,/-I/co(i-o) i a commonly 
used measure of star formation efficiency. 

Trend with nuclear separation 

Shown in Figure 0] is a plot of log(r3i) vs. the logarithm 
of the nuclear separation (log(sep , kpc)), for which we see a 
weak anticorrelation (r = —0.52). We observe no correlation 



between log(Lna) and the logarithm of the nuclear separa- 
tion (Figure[5| ■ This lack of correlation is unsurprising, given 
that log(Lfiri) is expected to depend on a number of differ- 
ent factors (e.g. the total extent of the star forming regions, 
dynamical mass of overall merger), which in turn are not 
simply, if at all, related to the progression of the merger. In 
contrast, by looking at the CO(3-2)/CO(l-0) line ratio we 
are in effect independent of the overall mass of the merger 
system. If we assume, therefore that r3i is a measure of the 
extent of excitation of the molecular gas present in the merg- 
ing galaxies, then the weak anti-correlation (r = —0.52) is 
consistent with the gas excitation increasing as the merger 
progresses. 



Correlations between line ratio, L fiR and star formation 
efficiency 

A common quantity used to determine the star formation ef- 
ficiency of a galaxy or galactic region, is log(LfiR/Z'CO(i-o))> 
where we assume here that the far-infrared luminosity is a 
result of emission from dust heated by UV-photons emitted 
from star formation rather than AGN. Shown in Figure [S] is 
a plot of log(r3i) vs. log(Lfnt/Lco(i-o))- A weak correlation 
(r = 0.38) is seen, suggesting that there is an increase in gas 
excitation in galaxies that show an increased SFE. This is 
not surprising, as one might expect to observe higher gas 
excitation in regions where star formation is proceeding vig- 
orously. No significant correlation is seen between log(r3i) 
and log(LfiR) (Figure[7l r = 0.21). 

Our observations and analysis hint at an increase in 
excitation, as traced by CO(3-2), with the temporal pro- 
gression of the merger. The uncertainties within our current 
data along with the relatively small number of sources per 
unit change in LfjR and nuclear separation make it difficult 
to state this with a sufficient degree of confidence, however. 
For the objects larger than a beamwidth, wide-angle jiggle- 
mappi ng with the newly commissioned HARP focal plane 
array (|Smith et al.i r2003) will enable the full extent of the 
CO (3-2) emission from the merging galaxies to be captured 
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0.5 1 1.5 

Log (separation (kpc)) 

Figure 4. log(r3i) vs. the logarithm of the nuclear separation 
for the detected objects in the LIG sample. Crosses indicate ob- 
jects observed with a single pointing, whilst open circles denote 
objects observed with multiple pointings. The error bars reflect 
1 — a errors in the CO(3-2) line intensities, and exclude ~ 20% 
absolute calibration errors. The best fit line, weighted by the ran- 
dom errors, is shown for illustration, and the Pearson correlation 
coefficient is -0.52. 



-A "+ 



Log( L Jm / L co(1 _ 0) ) 

Figure 6. A plot of log(Lfjf l /LQQ( 1 _Q)) vs. log(r3i) for detected 
objects in the LIG sample, with symbols defined as in Figure [4] 
The best line fit to the data, weighted by the random errors, is 
given and the Pearson correlation coefficient is 0.38. 



x O O 



0.5 1 1.5 

Log (separation (kpc)) 

Figure 5. A plot of log Lfijj vs. the logarithm of the nuclear sep- 
aration for the detected objects in the LIG sample, with symbols 
defined as in Figure[4] The Pearson correlation coefficient is 0.06. 



more accurately, enabling a more accurate line ratio to be 
determined. 



4.2 Contribution to the 850 /im continuum flux 

An 850 (jm / 450 continuum emission survey has recently 
been completed (|Mortier fc Is aak 2009) for the galaxies in 
the Gao and Solomon sample in order to map the the distri- 
bution of the dust emission and determine dust masses and 
temperatures. This has been mad e with SCUBA, a bo lome- 
ter array receiver on the JCMT (|Holland et al J 1 19991 ). The 
SCUBA 850 camera has a total bandwidth of ~30 GHz 
set by a bandpass filter with a centre frequency of 347 GHz. 
If the redshift of an object is such that the CO (3-2) line 
falls within the bandpass of this filter, the CO (3-2) line 
emissi on will contaminate the cont inuum emission from the 
dust (|Papadopoulos fc Alien! |2000h . The CO(3-2) observa- 
tions made here are thus very useful in determining the con- 



Log L fffi 

Figure 7. A plot of log(r3i) vs. log(Lfjj^) for the detected objects 
in the LIG sample, with symbols defined as in Figurc[4] The best 
line fit to the data, weighted by the random errors, is given and 
the Pearson correlation coefficient is 0.21. 



tribution of the CO (3-2) line emission to the 850 /im flux 
continuum flux measured using SCUBA. 

The CO (3-2) contribution to the measured SCUBA 
850 /im flux can be calculated using, 



S(V) 



2kv 2 L 



g(v) 



[Wm" 2 Hz _1 beam _1 ](3) 



ijSeaauist et al.l l2004h where vl is the line frequency, fit, 
is the main-beam solid angle and 7co(3-2) — / Jmb dV is 
the velocity integrated main-beam brightness temperature 
in Kkms -1 . For g(V), the bandpass profile of the filter 
expressed as a function of the recession velocity V = cz, 
we have adopted the same dual Gaussi an form used by the 
SLUGS survey in lSeaouist et all l|2004h . 

Table[3]shows the CO (3-2) contribution to the measured 
SCUBA 850 fim flux for the observed objects in the sample. 
As the SCUBA 850 fim fluxes are measured for the sample, 
these data will enable the correction for the CO (3-2) con- 
tamination of the dust emission to be made. These corrected 
850 /J,m fluxes will be published in a forthcoming paper - pre- 
liminary results show that the mean CO (3-2) contribution 
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Object 


CO(3-2) contribution 




(mjy) 




(Filter corrected) 


00057+4021 


2.7 


Arp 236 


91.4 


01077-1707 


14.9 


III Zw 35 


12.6 


Mrk 1027 


43.7 


02483+4302 


0.6 


UGC 2369 


13.8 


03359+1523 


7.4 


04232+1436 


0.0 


Arp 55 


21.8 


10039-3338 


9.6 


10190+1322 


0.0 


10565+2448 


9.9 


Arp 299 


103.1 


13001-2339 


27.2 


Arp 238 




NGC 5256 


27.1 


14348-1447 


0.0 


Arp 302 


37.4 


Mrk 848 


9.0 


Mrk 673 


10.7 


NGC 6090 


38.7 


17132+5313 


1.6 


NGC 6670 


30.6 


19297-0406 


().() 


20010-2352 


2.2 


II Zw 96 


15.6 


NGC 7592 


<4.1 



Table 3. Contributions of the CO(3-2) line to the 850 (im SCUBA 
continuum flux. 



for the sample is 22% and the median contribution is 17% 
ijMortierl 120091 . p riv. comm.). An 8 50 /im map of Arp 236 
has been made l|Fraver et al.l Il999h and the total flux at 
850 /im is reported at 273 mjy. The CO(3-2) contribution 
to this flux will be around 91.4 mjy (33%). These results 
demonstrate the importance of determining the correction 
both when fitting dust emission models to observed fluxes 
to determine dust temperatures, and particularly when de- 
termining total dust masses, whi ch are usually take n to be 
proportional to the 850 /xm flux l|Dunne et al.ll200oh . 



5 CONCLUSIONS AND FURTHER WORK 

We have observed the CO (3-2) emission from a sample of 
IR-selected merging galaxies using the single-pixel 350 GHz 
receiver at the JCMT. We have found that the CO(3-2) to 
CO(l-O) line ratio, rsi, is less than unity for all but one ob- 
ject in our merging sample with a sample average of r3i = 
0.47. LVG modelling indicates that this is consistent with 
T fe =(30-50) K, n — 10 3 cm- 3 and A = [CO/H 2 ]/(dV/dR) = 
10~ 5 (kms^pc- 1 )- 1 , although warmer (T k ~ 80-110 K) 
and more diffuse (n~3x 10 2 cm~ 3 ) conditions for the gas 
phase are not ruled out. We have calculated the CO(3-2) 
contribution to the 850 fim continuum flux for the galaxies 
in our sample. We have found that this contribution can be 
a large fraction of the 850 /im flux underlining the impor- 
tance of making CO(3-2)-based flux corrections when deriv- 
ing SED fits or dust masses from 850 /im fluxes. Our cor- 
rection will be presented in a forthcoming paper presenting 
the 850 /im SCUBA data for this galaxy set. 

The spread of line ratios observed in our merging sam- 



ple was broadly similar to that measured in both an IR- 
selected subsampl e of the SCUBA Local Universe Galaxy 
Survey (SLU GS) (lYao et al.l 120031). an d a sample of nor- 
mal galaxies (|Mauersberger et al.l l 1999) . We observe a weak 
anti-correlation between log(r3i) and log(sep; kpc) in our 
sample, a result consistent with gas excitation increasing 
as the merger progresses. A plot of line ratio versus star 
formation efficiency (log(r3i) vs. log(I/fiR,/L C o(i-o))) also 
shows a weak correlation, as gas excitation increases with 
increased star formation efficiencies. These correlations are 
suggestive of an increase in the excitation of the molecular 
gas, as measured by the CO(3-2)/CO(l-0) line ratio, with 
the progression of galactic mergers. However, studies of line 
ratios of these kind are often unavoidably affected by uncer- 
tainties in line luminosities arising both from random errors 
and systematic errors in the absolute flux calibration, which 
can be as high as 20%. Uncertainties also arise from assump- 
tions that often must be made about the angular size and 
shape of the CO (3-2) emission regions relative to the tele- 
scope's beam. These latter uncertainties will hopefully be 
reduced as more small sources {Da ~ a few arcseconds) are 
mapped interferometrically, and as more extended sources 
{Da > 10") are mapped with the next generation of multi- 
element focal-plane heterodyne arrays, such as HARP-B on 
the JCMT. 

Our future work includes an LVG analysis in which we 
will combine our CO (3-2) data with HCN line measurements 
for the same galaxy set to systematically probe the gas prop- 
erties across the merging sequence. We also intend improve 
our total luminosity estimates for our more extended sources 
by mapping them with the newly commissioned HARP-B 
focal-plane array. 



APPENDIX A: SYSTEM ATICS AND ERROR 
ANALYSIS 

As well as the systematic calibration errors of around ±20%, 
there are random errors arising from noise across each spec- 
tral channel. This noise leads to errors in determining the 
integrated spectral line and in d etermining the overall sub- 
tracted baseline level. Following [Gaol (1996), these random 
errors were combined using 

°I = \/4L + <sc = aT lmB VAWAV/Vl - AW/W.(A1) 

where aT TUla is the RMS noise on a single velocity channel 
of width AV, AW is the full width to zero intensity of the 
line and W is the total velocity coverage of the observation. 
These random errors in the velocity integrated main-beam 
brightness temperature, ol, are presented separately from 
any overall calibration errors. 



APPENDIX B: DERIVING THE GLOBAL 
CO(3-2) LUMINOSITIES FROM MULTIPLE 
POINTINGS 

For galaxies with a projected nuclear separation ~ 14" the 
procedure followed for combining the two measured inten- 
sities was similar to that used by IZhu et al.l (|l999h . Each 
galaxy was modelled as a Gaussian brightness distribution 
of diameter Do and inclination a with a position angle 



10 J. Leech, K. G. Isaak, P. P. Papadopoulos, Y. Gao, G. R. Davis 



Table Bl. Details of the combination of the observations for the galaxies observed with two pointings. For each galaxy the model with a 
small emitting region (1") was chosen, and the value of /total marked with an asterix used in calculating the final CO(3-2) luminosity. 



Object 


Pointing 1 


i«l 


D ,l 


PAi 


Pointing 2 


e*2 


£>(,, 2 


PA 2 


h 


12 


1 01 


hyl 


^TOTAL 


Arp 55 


(0,0) 





1 





(+10, +7) 





1 





28.1 


18.3 


26.0 


14.2 


40.2 


Arp 55 


(0,0) 





7 





(+10, +7) 





7 





28.1 


18.3 


30.8 


15.5 


46.3 


Arp 55 (corrected) 


(0,0) 





1 





(+10.+7) 





1 





28.1 


18.3 


24.4 


20.7 


45.1* 


NGC 5256 


(0,0) 





1 





(-7,-9) 





1 





18.4 


29.0 


13.1 


26.5 


39.6* 


NGC 5256 


(0,0) 





7 





(-7,-9) 





7 





18.4 


29.0 


13.9 


31.6 


45.5 


II Zw 96 


(0,0) 





1 





(+8,-9) 





1 





13.7 


17.9 


11.0 


16.2 


27.2* 


II Zw 96 


(0,0) 





7 





(+8,-9) 





7 





13.7 


17.9 


12.3 


19.0 


31.3 



PA. This procedure allowed the intrinsic intensities for each 
pointing 7oi and 7o2 to be derived from the measured inten- 
sity for each pointing 7i and 72. For II Zw 96, NGC 5256 
and Arp 55 detections were observed in both pointings and 
the parameters used in the source models are shown in Ta- 
ble lBll Objects where a line was only detected in one of two 
pointings were treated identically to the objects for which 
a single pointing was made, with the assumption that the 
CO (3-2) emitting region is small. 

The optical images of Arp 55, NGC 5256 and II Zw 96 
(Figure [Cll (w).(q') and (t)) suggest collisions between two 
face-on galaxies, and the optical brightness distributions are 
near circular so values of a — and PA = were chosen for 
both components. 

For each object, the intrinsic CO(3-2) intensities 7oi, 
7o2 were calculated for two cases: (i) assuming that the half 
power width of the CO (3-2) emitting regions were both small 
compared to the 14" beam (Do = 1")> and (ii) assuming 
the half power width of the CO (3-2) emitting regions were 
around one half of a beam size (Do = 7"). This gave an 
assessment of the sensitivity of the model to the unknown 
beam filling factors. The results for the two differing source 
sizes are shown in Table [Bll For calculating the final com- 
bined CO (3-2) intensity, the first model (i) corresponding 
to two small emitting regions was chosen in keeping with 
the assumptions made for the objects observed with a sin- 
gle pointing. Note that for Arp 55, it was found that the 
galaxy was inadvertently observed with a map centre of 
09:15:54.9 +44:19:54.4 (J2000) rather than the required co- 
ordinate of 09:15:54.7 +44:19:51 (see Figure [CiTwD. We 
thus incorporated this additional offset of (+1.9", +3.4") in 
each of the two pointings into our model, leading to an es- 
timated total flux of 45.1 K km/s (see the "Arp 55 (cor- 
rected)" row in Table [Bll above). 



APPENDIX C: CO(3-2) SPECTRA AND 
POINTINGS 

The main beam efficiency-corrected spectra together with 
the pointing centres superimposed on a Digital Sky Survey 
optical image, are shown in Figure [Cll 
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Figure CI. CO (3-2) spectra and optical images (wit h pointings a nd beam size). The optical image wavelength is 645 nm (band 103aE) 
(taken from the STScI Digitized Sky Survey (DSS) (Kent 1994), http://archive.stsci.edu/cgi-bin/dss_form ). Where the object 
was observed with both narrow (900 MHz) and wide (1800 MHz) spectrometer bandwidths, both spectra are show, with the wideband 
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Figure CI - continued 
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Figure CI - continued 
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Figure CI - continued 
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Figure CI - continued 
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